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ABSTRACT 

The discovery of numerous non-nuclear X-ray point sources with luminosities L > 
Iq39 gj.g g-i jj-^ several starburst galaxies has stimulated speculation about their nature 
and origin. The strong variability seen in several sources points to massive black 
holes as the central engines. If the flux is isotropic, the luminosities range up to 
?a 10^^ erg s~^, implying masses of A/>10'^A/q if the luminosity is sub-Eddington. 
Here we explore a model for these sources. We suggest that in some tens of percent 
of globular clusters a very massive black hole, M >50 Mq, is formed. This black hole 
sinks in < 10^ yr to the center of the cluster, where in the ^ 10^° yr lifetime of the 
cluster it accretes ^ 10'^ M©, primarily in the form of lighter black holes. Unlike less 
massive black holes in binaries, which are flung from clusters by recoil before they can 
merge gravitationally, a > 50 Mq black hole has enough inertia that it remains bound 
to the cluster. We suggest that '^10'^ Mq black holes may be common in the centers 
of dense globular clusters, and may therefore exist in some tens of percent of current 
globulars. If the cluster later merges with its host galaxy, accretion from young star 
clusters in molecular clouds by the black hole can generate luminosity consistent with 
that observed. We also consider the detectability of massive black holes in globular 
clusters with gravitational wave detectors such as LISA and LIGO, and speculate on 
future observations that may test our predictions. 

Key words: accretion, accretion disks - binaries: close - black hole physics - galaxies: 
starburst - globular clusters: general. 



1 INTRODUCTION 

In the last few years evidence has mounted for a popu- 
lation of spatially unresolved sources in starburst galaxies 
that have fluxes corresponding to isotropic X-ray luminosi- 
ties in excess of ~ 10'^^ erg s"^ (e.g., Fabbiano 1989; Fab- 
biano, Schweizer, & Mackie 1997; Colbert & Mushotzky 
1999; Zezas, Georgantopoulos, & Ward 1999; Kaaret et al. 
2001; Matsumoto et al. 2001; Fabbiano, Zezas, & Murray 
2001). The variability of these sources on time scales as short 
as 10* s (Matsumoto et al. 2001) suggests that they are 
accreting compact objects, probably black holes. If beam- 
ing of emission is moderate to negligible, the brightest of 
these objects (with L « 10*^ erg s~^) must have masses 
A/ > 10"^ A/0 so that the luminosity is below the Edding- 
ton luminosity Le = 1.3 x W^^ (M/IMq) erg s"^ Their 
non-nuclear locations imply an upper mass limit as low as 
10^~® Mq, otherwise dynamical friction would have caused 
them to sink to the center of their host galaxies (Kaaret 
et al. 2001). The inferred masses are much larger than the 
~ 10 — 20 Mq masses of black holes usually thought to arise 



from the evolution of present-day stars (Fryer 1999) , yet are 
much smaller than the ~ 10®~^ Mq masses of black holes 
that exist in the centers of many galaxies (e.g., Merritt & 
Ferrarese 2001). The origin of such objects is an intriguing 
mystery. It has been suggested that the emission is beamed 
(King et al. 2001), so that the accreting black holes have 
ordinary stellar masses, or that these black holes formed di- 
rectly via evolution of Population III stars (Madau & Rees 
2001). 

Here we propose another scenario, in which a black hole 
grows significantly by accretion. We show that for black 
holes with initial mass M < 10^ Mq, accretion from the ISM 
is too slow to increase the mass significantly. Accretion from 
stars in a dense cluster is more promising (see also Taniguchi 
et al. 2000), but the cluster must last several billion years for 
significant growth, which points to globular clusters. How- 
ever, a ~ 10 Mq black hole sinks to the cluster core within 
< 10^ yr, where it forms a binary or exchanges into one 
via three-body processes. Further three-body interactions 
tighten the binary, but each tightening causes the binary 
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to recoil. Before the binary becomes hard enough to merge 
rapidly via gravitational radiation, the recoil kicks it com- 
pletely out of the cluster (Sigurdsson & Hernquist 1993). 
Thus, a light black hole usually cannot increase its mass by 
a large factor despite being in a high-density stellar environ- 
ment. 

In contrast, as we show in § 2, a ^50 Mq black hole is so 
massive that recoil does not eject it from the cluster before 
the binary merges. Hence, the process can be repeated and 
in a dense cluster a > 50 Mq black hole can increase its mass 
to ~ 10^ Mq. This has some similarities to one of the ways 
in which supermassive black holes may grow in the nuclei of 
galaxies (see, e.g.. Prank & Rees 1976; Sigurdsson & Rees 
1997). 

We suggest that ~ 10^ Mq black holes are generated in 
this fashion in some tens of percent of clusters. While in the 
clusters, accretion onto the black hole is usually electromag- 
netically quiet, since most accretion is of compact objects 
that are swallowed whole. However, if the cluster merges 
with the disk of the host galaxy, as may occur for roughly 
half of globulars (e.g., Gnedin & Ostriker 1997; Gnedin, Lee, 
& Ostriker 1999; Takahashi & Portegies Zwart 1998, 2000), 
then after the globular has settled into the moan motion of 
the disk and has been disrupted by galactic tidal fields, the 
massive black hole is released. If it later enters a molecular 
cloud, the extra mass of the hole may speed up gravitar 
tional contraction and star formation, resulting in the black 
hole being at the center of (and accreting from) a cluster of 
stars. Por starburst galaxies, where the number of molecular 
clouds and volume-averaged density of interstellar gas may 
be a few times the corresponding values in normal galaxies, 
several such sources can be detectable at any given time. 

In § 2 wo estimate the growth rate of such black holes 
and the duty cycle of high-luminosity accretion (and hence 
detectability in globular clusters). In § 3 we calculate the 
expected accretion rate once in the disk of the host galaxy. 
In § 4 we discuss the implications of this model for the up- 
coming generation of gravitational wave detectors, and we 
present our conclusions and future observational tests in § 5. 

2 FORMATION AND GROWTH OF BLACK 
HOLES IN GLOBULAR CLUSTERS 

2.1 Birth of the initial >5OM0 black hole. 

Substantial growth of a black hole by mergers in a glob- 
ular cluster requires that it have a high enough mass, 
M > 50 Mq , to avoid being kicked out by recoil during three- 
body interactions (see § 2.2). How can such black holes form? 
We identify three possible mechanisms. The first is direct 
formation from a massive star. A number of recent numeri- 
cal studies indicate that if a star has a mass M > 30 — 40 Mq 
when its core collapses then the resulting explosion may 
not have enough energy to drive off the outer layers of 
the star (e.g., Pryer 1999). If so, and if the star has low 
enough angular momentum, then direct fallback from the 
star may produce a black hole with a mass close to that 
of the original star (Pryer 1999). Por a star with an ini- 
tial mass of ~ 50 — 100 A^q, the question is then whether 
stellar winds and mass ejection during the main sequence 
and giant phases are enough to reduce the mass well be- 
low the ~ 50 Mq threshold. The amount of mass lost in a 



wind is challenging to determine, but some recent calcula- 
tions suggest that, especially for stars in a low-metallicity 
environment such as a globular cluster, the wind loss is mi- 
nor; for example, from Table 3 of Vink, de Koter, & Lamers 
(2001), the mass loss rate for a 40 Mq star with metallicity 
Z = ^Zq is less than 1O~^M0 yr~^, implying a total of 
less than 1 Mq during the lifetime of the star. The number 
of such stars is difficult to estimate, given the poor statis- 
tics. Observations of stars in young clusters imply that for 
M > IMq the Salpeter mass function dN/dM oc M'"^'^^ 
is a decent approximation, with no hint of a breakdown for 
stars even as massive as ~ 100 M© (Massey 1998). This 
would imply ~ 10^* of stars have M > ^QMq, and hence 
that a globular cluster will potentially have tens of this mass. 

The second candidate mechanism is a variant of the 
first. In a dense stellar environment the most massive stars 
sink to the center of a cluster. If this happens more rapidly 
than the stars evolve off the main sequence (e.g., a few mil- 
lion years for a 30 Mq star), the stars may then undergo fre- 
quent collisions and mergers (Portegies Zwart et al. 1999), 
which produce more massive stars. Recent analyses (Pryer 
1999) suggest that when these massive stars undergo core 
collapse, they may leave 50— 100 Mq, black holes, even when 
mass loss from strong stellar winds is included. 

The third possibility is that if hundreds to thousands of 
~ 10 Mq black holes are formed in a cluster, there will be 
enough interactions that a few will merge with each other, 
forming a (smaller) set of ~ 20 Mq black holes, and so on, 
until the ~ 50 Mq threshold is reached. In this scenario, the 
overwhelming majority of small black holes are ejected from 
the cluster; it is only the fortunate minority that can merge. 
This mechanism has apparently not been explored in detail. 

2.2 Growth of >5OM0 black holes in cluster 
cores. 

Black holes formed from supernovae are much more massive 
than the average mass (m) w 0.4 — 0.6 Mq of stars in glob- 
ular clusters, hence they sink rapidly to the center of the 
cluster. A black hole of mass M sinks to the center in a time 
~ (jn)/M times the core half-mass relaxation time tr- Por 
typical globular clusters, tr ~ 10^~^ yr. Thus, a > 50 Mq 
black hole will sink to the center within < 10^ yr, shorter 
than the other time scales of interest. 

The rate of encounters in the core obviously depends on 
the number density of field black holes in the vicinity. Nu- 
merical simulations of globular clusters (Sigurdsson & Phin- 
ney 1995) show that in the cores of clusters in thermal equi- 
librium with number densities n > 10'^ pc~^ encounters of 
binaries with single stars are dominated by the most massive 
field stars, in this case black holes. This is because in ther- 
mal equilibrium the massive stars have a low velocity, hence 
they sink to the center of the cluster and have a high num- 
ber density oven if their total numbers arc comparatively 
small. Por example, in a multimass King model the scale 
height of a stellar species of mass m scales as ((m)/m)^''^ 
(see, e.g., Sigurdsson & Phinney 1995). Thus, black holes 
with typical mass 10Mq/0.4Mq = 25 times the average 
mass of a main sequence star have a scale height 1 /5 of the 
scale height of the typical star, and hence a number density 
enhanced by a factor of ~ 100. Therefore, if their number 
fraction averaged over the cluster is more than ~ 1%, their 
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core number density is comparable to or greater than that of 
visible stars. In addition, black holes may settle gravitation- 
ally with enough rapidity that they are no longer in thermal 
equilibrium with the lower-mass stars (e.g., Meylan 2000); 
in this case, their number density in the core and hence all 
encounter rates are correspondingly increased. For simplic- 
ity we assume that the core number density of black holes 
equals that of main sequence stars, but there is considerable 
uncertainty in this number. 

Wc focus on clusters that arc at or near core collapse, 
and adopt lO^ne pc~^ « 3 x lO'^'^ne cm~^ as a fiducial 
number density of normal black holes near the massive black 
hole. Of the 56 Galactic globular clusters studied by Pryor 
& Mcylan 1993 (out of a total of 143; sec Djorgovski & 
Meylan 1993), 6 had central mass densities greater than 
10^ Mq pc~^ and 22 had central mass densities greater than 
10^ Mq pc~^. The average stellar mass is less than 1 Mq 
(Pryor & Meylan 1993), so the number densities arc corre- 
spondingly higher. We show in § 2.2.2 that if na >0.1 the 
cores of clusters are favourable for the growth of black holes 
to masses ~ 10^ Mq, so this implies that some 40% of Galac- 
tic globulars, or « 60, could harbour such black holes. Note 
that this could be a substantial underestimate; for our pur- 
poses the important quantities are the masses and densities 
of clusters at birth, but to be conservative we are guided by 
the current properties of Galactic globulars. As mentioned 
by Portegies Zwart & McMillan (2000), zero age globulars 
may have had ~3 times the mass and ~ 1/3 of the virial 
radius of current globulars, and hence processes leading to 
the formation of massive black holes are likely to have been 
significantly more efiacient in the past. 

While in the high stellar density environment of the 
core, the black hole undergoes a variety of interactions that 
increase its mass via merger or accretion. If there is a suffi- 
cient fraction of binaries in the core (quantified in § 2.2.1), 
the massive black hole exchanges into a binary and ulti- 
mately merges. If the number of binaries is low, the black 
hole can capture a field star through tidal interactions (for 
a main sequence star) or by gravitational radiation (for a 
compact object). Three-body interactions strongly favour 
the resulting binary being composed of the two most mas- 
sive stars of the three that interact. The most massive field 
stars could be ~ 10 Mq main sequence stars within the first 
~ 10® yr of the cluster or, later on, could be ~ 10 M© black 
holes (if these are abundant), neutron stars, or white dwarfs. 
Further interactions tend to harden the binary, but they also 
tend to eject the third, field star. If instead accretion is dom- 
inated by captures of field compact objects by emission of 
gravitational waves, the compact objects are accreted with- 
out ejection. We now quantify these scenarios, starting with 
the demarcation between binary rich and binary poor. 

2.2.1 Division between binary-rich and binary-poor 
duster cores 

Which occurs more often, direct capture or binary exchange? 
A hard binary is defined to have an orbital binding energy 
greater than the kinetic energy of a typical field star. As- 
suming a velocity dispersion of 10'''t)ms,(i cm s^"*^ for 0.4 Af© 
field stars, this implies an orbital separation of ~ 10^^ cm for 
black hole binaries of total mass ~ 10 M©. Thus, if a mas- 
sive black hole encounters such a binary with a periastron 



rp less than 10^^ cm, there is likely to be a strong interac- 
tion with am exchange. This is w 10* times the 10^^ cm 
periastron necessary for direct capture (see § 2.2.3); since 
in the strong focusing limit the cross section scales linearly 
with the periastron, (TcoII ~ TTrp{2GM/v^) (where Vao is 
the velocity of the field star at infinity), this implies that if 
the binary fraction is much larger than ~ 10~* in the core, 
a ~ 50 Mq black hole will interact with binaries far more 
frequently than it will capture single stars. 



2.2.2 Binary-rich clusters 

Now consider the case in which core binaries are plentiful by 
the above definition. This may be true even if the primordial 

binary fraction is low, due to processes such as three-body 
binary formation (e.g., Lee 1995). If a hard binary is formed, 
subsequent close three-body interactions pair the massive 
black hole with the highest-mass species present in abun- 
dance (a consequence of exchange interactions), and tend 
to harden the binary further (Hcggie 1975). However, for a 
stellar black hole of "normal" mass M < 10 Mq, the recoil 
of the binary that accompanies its hardening ejects the bi- 
nary from the cluster before the binary separation becomes 
small enough that merger by gravitational radiation occurs 
faster than the next binary encounter (Sigurdsson & Hcrn- 
quist 1993). This may not be true for a more massive black 
hole. 

Ejection or retention of binaries. — To estimate the 
threshold mass, we adopt the treatment of Portegies Zwart 
& McMillan (2000). We are interested here in binaries with 
very unequal mass, but we assume for simplicity that the 
field black hole has a mass comparable to the mass of the 
lighter hole in the binary. For three equal-mass stars the 
binding energy is usually increased by ~ 20% by a strong 
interaction (e.g., Heggie 1975; Sigurdsson & Phirmey 1993), 
but for high mass ratio binaries Quinlan (1996) finds that 
in a given hardening interaction the binding energy of the 
binary typically increases by a fraction ~ 0.2(m/A/) = 
0.04mioMgQ^, where the mass of the large black hole is 
50Af5o Mq and of the smaller black holes is lOmio Mq. The 
fraction of the change in binding energy carried away by bi- 
nary recoil also depends on the mass ratio. For three equal- 
mass objects, conservation of momentum implies that the 
kinetic energy of recoil of the binary is approximately a third 
of the change in binding energy. If the binary is much more 
massive than the third star, M ^ m, then the kinetic energy 
of the binary is instead roughly m/M times the change in 
binding energy. Let Dose be the escape velocity from the core, 
where typically Vesc ~ 50 km s~^ for a dense cluster (Web- 
bink 1985). Then the binary binding energy that typically 
produces recoil at the escape velocity is 

Sb.n,i„ « 5{M/mf (^Mve'sc) « 2 X W^°Mlom^o^ erg . (1) 

This needs to be compared with the binding energy such 

that gravitational radiation causes a merger before the next 
three-body interaction. The merger time for two stars with 
total mass M and reduced mass /i in an orbit of semimajor 
axis o and eccentricity e is (Peter 1964) 

Wge = 3 X lO®M|(MM')-i(a/i?0)*(l - ey^^yi. (2) 
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From numerical simulations of three equal-mass black 
holes (Portegies Zwart & McMillan 2000), the typical eccen- 
tricity distribution of binaries after the interaction of three 
equal-mass objects is P(e) = 2e (with a slight excess at 
high eccentricities). With this distribution, the average ec- 
centricity weighted over merger time is e « 0.7. We note, 
however, that the strong dependence of merger time on ec- 
centricity means that deviations from this distribution can 
have large effects on the rate of growth of large black holes in 
globular clusters, the detectability of gravitational radiation 
from these systems, and other important quantities. Quin- 
lan (1996) finds that when perturbations are small (as they 
are when the mass ratio is high) the eccentricity of a hard 
binary is increased steadily by the hardening process un- 
til gravitational radiation becomes important. If so, merger 
times are decreased and conditions are more favourable for 
the growth of ~ 10^ Mq black holes in the centers of globu- 
lar clusters. However, in the presence of finite perturbations 
the outcome is not as clear. An important future project is 
therefore a numerical study of the growth of > 50 Mq black 
holes in globulars, with a full treatment of the evolution of 
the eccentricity. In the meantime, we assume conservatively 
that -P(e) = 2e regardless of the amount of hardening that 
has occurred, until the point that gravitational radiation is 
significant. 

The typical semimajor axis for a merger time Tmcrgc ~ 
lO^Te yr and M ^ m is then 



o inll l/4»#l/2 1/4 

a ^ 3 X 10 Tg ' Mgp m^Q cm 



(3) 



The binding energy is 

Bbind.mergc = GMm/a 4 X lO^Vg" Afgf 771?^'' crg . (4) 

Therefore, the ratio of the binding energy necessary to eject 
the binary from the cluster to the binding energy necessary 
for a fast merger is 



,min / -Ebind , merge 



11/4 1/4 
10 '6 



(5) 



If a:: > 1, mergers happen before ejection, so that the hole 
gains mass and remains near the mass source. If a; < 1, 
ejection happens first, so that although the binary might 
merge in less than a Hubble time it will do so well away from 
the high stellar density in the cluster, meaning that the black 
hole mass growth is stopped. This calculation shows that if 
Af > 50 Mq the hole will stay near the center of the potential 
and grow rapidly, whereas for much smaller masses the hole 
will be ejected quickly without additional mass growth. If 
the eccentricity increases as the orbit shrinks, the critical 
semimajor axis increases, and hence the ratio x increases. 
Therefore, eccentricity growth allows lower-mass black holes 
to grow by mergers without being ejected. 

Time scale for growth by three-body interactions. — 
The interaction cross section for a massive black hole 
is dominated by gravitational focusing, so that CTcoII ~ 
7rrp(2GAf/«^). The relevant pericenter distance depends on 
the type of encounter and the nature of the secondary ob- 
ject. For example, a 10 Mq BH companion with a typical 
eccentricity e — 0.7 and a semimajor axis ~ 5 x lO'^^Tg^'* cm 
will merge via gravitational radiation within lO^re yr. A 
third black hole passing within ~twice that distance will 
produce a strong interaction, so rp ~ 10^^ cm. In thermal 
equilibrium the velocity dispersion of stars is inversely pro- 
portional to the square root of their mass, so if the average 



mass of a main sequence star is 0.4 Af© and its velocity is 
VuiB = 10''ums,6 cm s~^ then black holes have a velocity 
dispersion v ~ ^nij^Q^^VmB — 2 x lO^m 
Then ctcoii « 10^°(rp/10^^ cm)Af5om 
time between interactions is Tint ~ l/(naVoo) ~ 2 x 



-1/2 



«nis,6 cm S . 

m^ , so the 



lO^'vms.en^^r'l^M^g^m'f^^^ yr, where rp = 10^^rp,i2 cm 
and the average in angle brackets is taken over the veloc- 
ity distribution of black holes (assumed to be Maxwellian). 
The binding energy scales with 1/a, as does the time be- 
tween interactions. Therefore, the orbit starts shrinking 
rapidly, but slows down. The total time required to reach 
a = 5 X 10^^ cm from an initially much larger semimajor axis 
is ~ 5(A//m)rint if a fraction ~ 0.2{m/M) of the binding en- 
ergy is removed in each close encounter. The orbital separa- 
tion for a fixed merger time scales as M^^^, so Tint ~ -^50^''^ 
and a ~ 50 Mq black hole will increase its mass by ~ 10 Mq 
in a total time ~ 
ing, after time t the 

t/ (s X 10** n;. 



5 X lO^nr^AT. 



-1/2 -5/2 

'10 



-5/2 

ho yr 



yr. Integrat- 

mass of the black hole is Afso ~ 
2 

. Within a Hubble time there 



is therefore time for substantial growth in the mass of large 
black holes in the centers of clusters that are at least moder- 
ately dense, with ng ^0.1. If, as argued by Quinlan (1996), 
three-body effects increase the eccentricity as the binary 
hardens, the rate of increase of black hole mass is enhanced 
significantly. 



2.2.3 Binary-poor clusters 

Now consider what happens when binary interactions and 
exchanges can be neglected. The effective capture cross sec- 
tion of a compact object of mass m by a large black hole of 
mass Af 3> m is (Quinlan & Shapiro 1989) 



85^ V GW^^Afi^A 



(6) 



Here Voo is the relative velocity at infinity and c is the speed 
of light. Numerically, cr « 2 x IQ^^'ml^J mII'\~^^^'^ crn^ , 
where Voa ~ lO^we cm s~^. The strong dependence on rel- 
ative velocity means that for equal number densities the 
slowest moving population will dominate the encounter rate. 
In thermal equilibrium, this means the black holes. Av- 
eraging the encounter rate {navoo) over the velocity dis- 
tribution (assumed Maxwellian), this implies ct « 2 x 
lO^^mJJ^^AfgQ^^Wj^^e'''^ cm^ (recall that Vao oc m'^^'^Vms). 
The periastron is only ~ 10^^ cm for such an encounter, im- 
plying a merger time of ~ 10** yr once the orbit is circularised 
(see equation (|^; circularisation also takes place rapidly). 
Thus, if a compact object is captured in this manner it is 
swallowed almost immediately, without further interactions. 
The encounter rate for a massive black hole is then 



/ \ o ir\ — -Lb II i n /f-L^/ ' 

nc = {navao) ~ 2 x 10 nenij^g Afgg' v. 



11/7 »,rl2/7^ -11/7^-1 



(7) 



The characteristic time scale to increase the black hole's 



mass by of order itself is then tacc ~ {M/m)u^^^^, or 
t ~ 2 X 10%-^m"^**/''M"^''''7)"'"'vr 



(8) 



Thus, for a dense cluster (ng ~ 1), or an especially high- 
mass black hole {M50 > 5) , significant growth can occur in a 
Hubble time even in the absence of binaries in the core. Note 
that the growth rate accelerates with increasing M, so once 
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growth is significant the massive black hole will consume the 
remaining field black holes quickly. The ratio of three-body 
to two-body encounter rates is 



2 X 



(9) 



Wo sec that unless there arc virtually no binaries in the core 
into which a massive black hole can exchange, such a black 
hole will encounter binaries far more often than it encoun- 
ters single stars. However, note that merger via three-body 
encounters may require tens to hundreds of interactions (see 
below), whereas merger via direct two-body capture is es- 
sentially immediate. Moreover, just as three-body rates are 
modified by gravitational radiation, the direct capture rate 
may be modified by deflections induced by three-body ef- 
fects. Numerical simulations are required to determine when 
the two-body growth rate exceeds the three-body growth 
rate of a massive black hole in a cluster, and how the two 
effects couple. 



2.3 Mass of central black holes in globular clusters 

If the escape velocity from the core is 50 km s-i (Web- 
bink 1985) and the field black hole is ejected with a ki- 
netic energy a factor w 0.2(m/M) of the binary binding 
energy, then field black holes will be ejected when the bi- 
nary orbital radius is about lO^^mio cm, independent of 
the mass of the larger black hole. At a fractional harden- 
ing of 0.2(m/M) per interaction, the number of interac- 
tions required to bring the radius down to w 3 x 10^^ cm 



In 



10^ 



/ln[l -F0.2m/M] 



^mio/(3 X 10"M5'/"m;/*) 

2bMsQrn~^^\n(mM~^^'^rn%'^) ^ 100 M50. Thereafter it 
merges quickly by gravitational radiation. Thus for each 
black hole accreted in this manner, w 100 Mbo are ejected, so 
the massive black hole can in principle consume ~ 0.1—1% of 
the smaller black holes that interact with it, ejecting the rest. 
Here again the precise eccentricity behaviour has a large im- 
pact; if high eccentricities are reached in the process of hard- 
ening, the binary merges when the semimajor axis is larger, 
and hence ejects fewer field black holes. The smaller black 
holes will also interact with and eject each other, further 
cutting down the mass supply. Portegies Zwart & McMillan 
(2000) find that interactions among the smaller black holes 
eject roughly half of them within ~ 2Gyr, and about 90% 
after several billion years. Thus, if a large black hole grows 
rapidly, it may be able to accrete a few tenths of a percent 
of the initial mass in black holes. How much mass would this 
be? 

If the initial mass function of stars in the cluster is the 

Salpeter IMF, dN/dM oc M"^'^^ above 1 Mq (and fiatter 
below one solar mass; see, e.g., Meyer et al. 2000), then the 
mass in stars above some threshold Mo > 1 Mq is 



M^dMocMo°-'' 
dM 



(10) 



If black holes form from stars with initial masses m > 
30 M0, this implies that ~20% of the initial mass might 
have been in stars that form black holes. Wind losses on the 
main sequence, mass loss during the supernova, and ejec- 
tion of some black holes by supernova recoil reduce this 
fraction. However, wind losses may be minimal for these 



low-metallicity stars (Vink, de Koter, & Lamers 2001), and 
supernova recoil of black holes is expected to be small com- 
pared to that for neutron stars (e.g., Sigurdsson & Hernquist 
1993), so mass losses could be unimportant. This is espe- 
cially true if stars of initial mass AI > 40 AIq produce failed 
supernova (Fryer 1999) and hence lose little mass. Given 
the initial total mass in black holes and the fraction lost to 
ejection, a large central black hole could accrete a fraction 
~ 10~^ of the initial mass of the cluster. Perhaps coinciden- 
tally, this is comparable to the ~ 0.5% ratio of supcrmassive 
black hole mass to bulge mass found for many galax;ies (e.g., 
Merritt & Ferrarese 2001). Then, an initial cluster mass of 
^ W^' Mr.), appropriate for the densest clusters, implies that 
the black hole could gain a few hundred to a thousand solar 
masses from just the smaller black holes. If there are multi- 
ple ~ 50 Mq black holes that can merge with each other, or 
if high eccentricities are reached, this number may become 
~ 10^ - 10^ A/c?. 

Once the small black holes are depleted significantly the 
main source of fuel becomes the next most massive stars in 
the core. If the black holes are consumed before > 1.5 Mq 
main sequence stars evolve, these will dominate the core. 
Otherwise, the next in line are neutron stars. These neutron 
stars have a number density comparable to or smaller than 
that of black holes (having a larger scale height and not be- 
ing overwhelmingly more numerous) and are also less mas- 
sive, so the characteristic time for mass gain becomes longer 
by a large factor, perhaps ~ 100 or more. Thus, some extra 
growth is possible, but not factors of many. It is therefore 
plausible that this mechanism would preferentially generate 
black holes with masses M ~ 10^ Mq. 

2.4 Detectability of black holes in clusters. 

Most globular clusters obviously do not have 10*^ erg s"'^ 
X-ray sources. If a substantial fraction contain high-mass 
black holes, why are they not emitting profusely? To answer 
this, consider what happens when stars of various types are 
accreted. Also, remember that the most massive type of star 
present in abundance will dominate the encounters. This is 
because binary exchange interactions tend to yield a binary 
consisting of the two most mgissive of the three stars in- 
volved, and a large number of such interactions are required 
to harden the binary sufficiently to produce a merger. If, for 
example, 100 binary interactions are required before merger, 
then the encounter rate with black holes must be < 1% of 
the total for the final binary to not include two black holes. 
Given that the lower velocities of black holes increase their 
interaction cross section, this would require severe depletion 
of the black holes in the core. If instead direct capture dom- 
inates, the higher mass and lower velocities of massive ob- 
jects give them a far higher cross section than the less mas- 
sive objects. Thus, mergers are primarily with black holes 
when they are present, thou uoutron stars, and finally white 
dwarfs and the remaining main sequence stars (which have 
M < 0.8 Mq for a typical current cluster age). 

A small black hole merging with a massive black hole 
will generate only gravitational radiation. The same is true 
with a neutron star. The critical orbital separation rtidc be- 
tween an object of mass M and one with mass m <C M and 
radius R, such that the less massive object is destroyed by 
tidal forces, is rude ^ 2.5{M/m)^^^R, or about 10^ cm for 
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M = 50 Mq and m — 1.5 Mq. Assuming that the massive 
black hole is not rotating significantly (plausible, since inter- 
actions are likely to be uniformly distributed in direction), 
the radius of the innermost stable circular orbit, inside of 
which an object will plunge without further loss of angular 
momentum, is i?isco = &GM/c^ = 5 x lO'^ M50 cm. The 
two radii are equal, rtide = -Risco, when M w 5Mq, and 
since -Risco/^'tide oc M^^^ more massive black holes swallow 
neutron stars even more easily. Therefore, a neutron star is 
not disrupted before merger. 

A white dwarf is disrupted, but the merger is very 
quick. A white dwarf with typical mass 0.6 Mq has a ra- 
dius R ~ 10^ cm (e.g.. Panel, Althaus, & Benvenuto 2000). 
The tidal disruption radius is then « IO^^Mjq^ cm, well 
outside the innermost stable circular orbit. The dwarf there- 
fore transfers mass as in a low-mass X-ray binary. However, 
at such separations the time to merge via gravitational ra- 
diation is short: roughly ten years for Mso = 1, scaling 
as M~^/^. Even if this time increases as the white dwarf 
loses mass, there is only a short time in which radiation can 
be detected. In addition, the initial mass transfer rate of 
~ 0.1 Mq yr"^ is 10* times the Eddington rate for a 50 M© 
black hole. At such high rates emission may be in the hy- 
percritical regime in which neutrino emission dominates over 
photon emission (Ifouck & Chevalier 1991; Brown, Lee, & 
Bethe 2000). Thus, white dwarfs, if accreted, will have a tiny 
duty cycle for observable emission, if indeed the emission is 
observable at all. 

A main sequence star of mass ~ 0.6 Mq has a radius 
~40 times larger than a white dwarf of the same mass. The 
separation at disruption is thus 40 times higher, and the 
gravitational merger time (40)* ~ 10^ times greater, than 
for a white dwarf. Thus, accretion of a main sequence star by 
a massive black hole will produce copious photon emission. 
At the Eddington rate a 0.6 Mq star will be consumed by 
a ~ 50 Mq black hole in ~ 2 x 10^ yr, smaller than but 
comparable to the hardening timoscalc. If a cluster black 
hole has exhausted its supply of compact objects with mass 
greater than the main-sequence turn-off mass, it is therefore 
possible that it may form an LMXB with a main sequence 
star. 

Accretion from gas in the cluster is unlikely to be de- 
tectable. Although much of the primordial gas in globular 
clusters may be swept out by supernovae, the winds from 
evolving red giants may be accreted. If we assume that at 
any given time 1% of the stars axe on the red giant branch, 
even a rich cluster with 10® stars will have only 10'' red gi- 
ants, which arc distributed within a typical half-mass radius 
of 10 pc. The mass loss rates for low-mass red giants are dif- 
ficult to determine but seem to be < 10~* Mq yr~' and may 
be significantly less for the metal-deficient stars in globular 
clusters (Dupree 1986). The winds come out in a large range 
of velocities broadly distributed around the escape velocity 
(Lamers & Cassinclli 1999), which is roughly 50 km s^^ 
at the 1 AU radius of a 0.8 Mq red giant. The winds take 
3x 10'^ cm/5 X 10® cm s~'=2x 10* yr to cross the volume, so 
that on average each star donates « 10~^ Mq in wind ma- 
terial, for a total of 10 Mq in ^(10pc)=' ^ 10*'^ cm^ The 
density is then 2 x 10"^'^ g cm^"^, consistent with the number 
density of ~0.1 cm~^ inferred from pulsar dispersion mear 
sures in 47 Tuc by Freire et al. (2001). Suppose the central 



black hole captures the gas via Bondi-Hoyle accretion, at 
the rate 

M « 2 x 10^^Mlop-2i [vl + c%\ g s-i , (11) 

where the hole has a velocity 10® ue cm s~^ relative to the 
cloud, and at infinity the density and sound speed of the gas 
are 10~^''p-24 g cm~^ and 10® ce cm s~', respectively. Then, 
the accretion rate is ~ 5 x 10^'^ g sT^ for a 10'' A/q black 
hole. This is, however, likely to be an overestimate by sev- 
eral orders of magnitude. The X-ray luminosity generated 
by accretion ionises and heats its surroundings (e.g., Mal- 
oncy, Hollenbach, & Tielens 1996) , so that even well beyond 
the Bondi-Hoyle capture radius the effective temperature is 
raised to T > 10® K by Compton heating. Pressure balance 
dictates that this heated gas must be less dense than its 
surroundings by a factor proportional to the temperature, 
so that the net accretion rate can be far lower than that 
indicated by a simple Bondi-Hoyle estimate. Accretion from 
any form of interstellar gas, molecular, atomic, or other- 
wise, is thus unlikely to produce the observed luminosity. It 
is therefore difficult to rule out the presence of a large black 
hole based on upper limits to X-ray emission. For example, 
Grindlay ct al. (2001) have recently observed the cluster 47 
Tuc with the Chandra X-ray Observatory. Even without in- 
cluding the possible quenching of accretion by preheating, 
they find that despite the lack of any detectable point source 
of X-rays within 3" of the cluster center, a black hole with 
a mass of up to 500 Mq could have escaped detection. 

A final way in which such a black hole could be de- 
tected in a globular would be through its long-range effect on 
other cluster stars. Current data arc ambiguous. For exam- 
ple, Sosin (1997) finds that in the post core collapse clusters 
M15 and MSO the radial density cusp of stars is more consis- 
tent with the presence of a ~ 10^ Mq black hole than with 
simple equipartition core-collapse models. However, there is 
only marginal consistency with the expected cusp in the ve- 
locity dispersion (e.g, Gebhardt et al. 1995). Observations 
of the pulsars in M15 suggest that there is no more than 
~ 700 Mq in black holes in the center of this core- collapsed 
cluster (Phinney 1993). Future observations will be required 
to place more stringent limits on the total black hole mass 
in various clusters. 



3 CAPTURE OF CLUSTERS BY THEIR HOST 
GALAXIES AND ACCRETION IN THE 
GALACTIC DISK. 

Globular clusters around a galaxy such as our own interact 
with the galaxy through tidal stripping and other events. 
Current estimates are that approximately half of the clus- 
ters eventually merge with the disk and eventually disperse 
(e.g., Gnedin & Ostriker 1997; Gnedin, Lee, & ©striker 1999; 
Takahashi & Portegies Zwart 1998, 2000). A massive black 
hole generated in the center of such a cluster is then released 
and can interact with the interstellar medium. An intrigu- 
ing possible effect of a massive black hole inside a molecular 
cloud is that it may help precipitate cloud collapse and star 
formation. Note that in a giant molecular cloud of density 
p = 10~^^ g cm~^, a 10^ Mq black hole contributes most 
of the mass out to a distance of a few parsecs. The natu- 
ral gravitational instability of the cloud would then be en- 
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hanced by the presence of the hole, and in a few free-fall 
times iff ~ l/VGp ^ 10''(p/10"^^g cm"^)"^/^ yr a star 
cluster could form. This is comparable to the crossing time 
of a giant molecular cloud by a massive black hole, assuming 
a relative velocity of a few km s~^ (appropriate if the black 
hole is nearly at rest relative to the local average, since gi- 
ant molecular clouds tend to have a random bulk velocity 
of 3-5 km s~^; see Dame, Hartmann, & Thaddeus 2001). As 
in a globular cluster, the most massive stars would sink to 
the center, where (as main-sequence stars) they would ac- 
crete onto the black hole and generate significant luminosity, 
in principle up to the Eddington luminosity. This may ac- 
count for the observation by Matsushita et al. (2000) that 
the brightest point X-ray source in M82 appears to be in a 
star cluster within a superbubble. 

Note, however, that a black hole with a roughly stellar 
mass M ~ 10 Mq cannot easily grow to the required ~ 
10^ M0 masses in such an environment. This is because the 
e-folding time for mass increase at even the Eddington rate 
is 6 X 10^ yr, comparable to or longer than the timescale 
on which young star clusters dissipate (Portegies Zwart et 
al. 2001), so that although a black hole may be luminous 
as it feeds off of massive stars, the amount of mass it can 
accumulate is relatively small. These black holes must grow 
elsewhere. 

How likely is it that a ~ 10^ Mq black hole is in a 
dense interstellar cloud at any given time, and what is the 
resulting prediction for the number of active sources at any 
given time? Consider first our Galaxy. It currently has ~ 150 
globular clusters, of which ~40%, or 60, have core densities 
ne > 0.1 (Pryor & Meylan 1993) and therefore produced 
massive black holes in their cores according to our model. 
We assume that a similar number have merged with the 
Galaxy in its lifetime. Molecular clouds with densities > 
lO'^ cm~^ occupy a volume fraction < 10^^. It is therefore 
not surprising that no massive black holes are currently in 
dense clouds in our Galaxy. 

In contrast, galaxies undergoing active star formation 
are thought to have surface densities of molecular gas that 
are higher than ours by a significant factor, perhaps ~ 10 or 
so (Taniguchi & Ohyama 1998; Gao et al. 2001). If the molec- 
ular clouds in star-forming galaxies are similar to those in 
our Galaxy, the volume fraction of dense clouds is therefore 
likely to be up by a similar factor, to ~ 10~^. That sug- 
gests that such a galaxy will have a few massive black holes 
currently in star clusters within molecular clouds, emitting 
> 10'^^ erg s~^. This is consistent with the number observed. 



4 DETECTION OF GRAVITATIONAL WAVES 
FROM A 10^ Mq black HOLE IN A 
GLOBULAR CLUSTER 

A major prediction of our model is that tens of percent of 
globulars have ~ 10^ Mq black holes. As discussed in § 2.4 
and § 3 these holes are usually electromagnetically quiet. 
However, a massive black hole in a globular cluster around 
our Galaxy is a potentially strong source of persistent grav- 



itational waves in the 10 



10 Hz sensitivity range of 



planned space-based interferometers such as LISA. In the 
final stages of inspiral some of these sources may be de- 
tectable with ground-based detectors such as LIGO II up to 



~ 1 Gpc away. The strength, detectability, and number of 
such sources depends on many uncertain parameters, such 
as the eccentricity distribution and the current number den- 
sity of lighter black holes in globulars. In this section we 
make rough estimates of the type and strength of signals 
that may be found. 

The dimensionless gravitational wave strain amplitude 
measured a distance r from a circular binary of masses M 
and m with a binary orbital frequency /bin is (Schutz 1997) 



h = 2^/^(4^)^/3^/,^i(fMm(M + m) 



-1/3 



1 



(12) 



Note that the frequency of the gravitational waves is twice 
the binary frequency. A black hole binary will emit grav- 
itational waves at all separations, but here we focus on 
the portion of the inspiral during which gravitational ra- 
diation has a significant infiuence on the orbital evolution. 

years, for a duty cycle of 



From §2, this may happen in roughly the last 10^ Mg, 
of a total of 5 X 10^ M"^^^ 



~ 20%. For a inspiral time of 10^ M^^^^'^ yr and a cen- 
tral density of lO^ne pc~"^, the semimajor axis is a ~ 



4 X 10 cmn,i 



1/4,^3/8 1/4 
' M, ' 



50 



SO that the binary orbital fre- 
quency is /bin = 4x lO'^Hz rig'^^M^^^'^^m^Q^^*. Numerically, 
the gravitational wave amplitude is therefore 

/i « 2 X 10"^^ne/*M5/Vi^''(10kpc/r) . (13) 

The sen sitivity of LISA (from the LI SA 



homepage, http://lisa.jpl.nasa.gov/science/science.html) in 
the range 10^* — 10^'^ Hz at a gravitational wave fre- 
quency /gw = 2/bin is approximately /imin = 3 x 

10-21 (/gw/10-* Hz)-^/2 ^ 5^^Q-22(^^.^^/^g-4jj^-)-5/2 f^j. 

a signal to noise ratio of 5 in an integration time of one year. 
Thus, the signal to noise ratio for most of the inspiral after 
the last three-body interaction is 

-3/16 lOkpc 



S/N = 5 h/hn 



5n, 



19/16^15/32 



6 



50 



(14) 

Therefore, black holes with masses greater than ~ 100 Mq 
can be detected with high significance in a year in a dense 
cluster within 10 kpc, and in a few years black holes 
with masses greater than ~ 10'^ Mq can be detected if 
ng > 0.1. The steeply improving sensitivity of LISA in the 
10"'' - 10"^ Hz frequency range means that gravitational 
waves from a neutron star or white dwarf being accreted are 
slightly easier to detect than black holes being accreted. If 
20% of mergers are in a favourable phase of the evolution, 
and if ^40% of Galactic clusters have central number densi- 
ties ne > 0.1 (Pryor & Meylan 1993), then ~ 10 cluster black 
holes will be in this phase and detectable if M > 10^ Mq . 

The above applies to the amplitude of waves emitted 
during most of the spindown. In addition, some fraction 
of sources will be in the final stage of their spindown, and 
hence will emit waves of higher amplitude. For a given pair 
of masses the amplitude scales as /^/f and the sensitivity 
threshold of LISA scales as f^^^^'^ , so the signal to noise 
scales as /bf/^ up to the « 10"^ Hz sensitivity peak of LISA. 
The fraction of sources with a given orbital frequency scales 
as the inspiral time, which is ~ a* ~ fhil^^ ■ Intermediate 
mass black hole binaries may therefore be detectable at large 
distances; for example, at the 15 Mpc distance of the Virgo 
cluster a fraction ~ 10"^ of such binaries are likely to be 
observed, for a total of perhaps tens of sources among the 
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~ 10^ galaxies in the cluster. Binaries in the last stage of 
inspiral will of course be visible to even larger distances. 
For example, from equation a 10 Mq black hole in a 
circular orbit around a 10"^ Mq black hole at a binary or- 
bital frequency of 5 x 10~^ Hz will spiral in within « 10 yr. 
The dimensionless strain amplitude at a distance of 10 kpc 
would be ft = 4 X 10"^^, compared to the LISA S/N = 5 
sensitivity oi h = 10"^"^ at 10~^ Hz. Such a signal would 
therefore yield a 1 year signal to noise of 10 out to 200 Mpc 
with LISA. Portegies Zwart & McMillan (2000) estimate the 
number density of globulars in the universe as 8Ah^ Mpc"^; 
if ~ 10% of these have 10'^ Mq black holes and core number 
densities ng « 1, we would expect ~ 10 such binaries to be 
detectable at any given time. 

The merger itself is at frequencies above the range of 
LISA. However, if an intermediate mass black hole in a 
binary within a few years of merger is found with LISA, 
the final inspiral could be anticipated and detected with 
ground-based detectors such as LIGO, because the fre- 
quency, phase, and location of the source would be known 
in advance. Unfortunately, the frequency of this final merger 
[few ~ 40 (A//10^ A/q)"^ Hz at the last stable orbit of a 
nonrotating black hole] is only marginally observable with 
LIGO-I, which has a lower frequency cutoff ~ 40 Hz (see, 
e.g., Brady & Creighton 2000). However, LIGO-II is ex- 
pected to be able to detect 10 Hz sources at 99% confidence 
for a 10 second integration and a dimensionless strain am- 
phtude oi h — 2 X 10"'^'^ (scaling from Figure 3 of Brady 
& Creighton 2000). If black holes at the low end of the 
masses we consider (100 Mq instead of 10^ Mq) are present 
in a significant number of globulars, their final merger is in 
the LIGO-II sensitivity band. From equation (^ , a 10 Mq 
black hole spiraling into a 100 Mq black hole spends ap- 
proximately 10 seconds at a gravitational wave frequency of 
10 Hz or higher. Using equation ( p^ ) as a rough guide, this 
signal would potentially be detectable out to a gigaparsec or 
more. Thus, advanced ground-based detectors may pick up 
gravitational waves from a large number of ~ 100 Mq black 
hole mergers in globular clusters. 



5 DISCUSSION AND CONCLUSIONS 

We suggest that a population of ~ lO'^ Mq black holes is gen- 
erated in globular clusters and then released into the disks of 
their host galaxies when the clusters are assimilated. Subse- 
quent accretion from young star clusters formed from molec- 
ular clouds then produces the luminosity observed from non- 
nuclear point sources in other galaxies. Such black holes, ac- 
creting at less than the Eddington limit, can produce the 
observed flux without strong beaming. 

King et al. (2001) have proposed instead that the black 
holes are of ordinary stellar mass but are strongly beamed. 
They argue that larger black holes would have difficulty in 
producing the required high-energy X-rays if the emission is 
from blackbody annuli in the disk, because the characteris- 
tic blackbody temperature for an Eddington accretor near 
the inner edge of the disk is kT ^ 2(M/M0)"^/* keV. How- 
ever, we point out that if a hot Compton corona exists near 
the black hole, as is usually invoked to explain the high- 
est energy emission, the temperature scales as M/r and is 
thus independent of mass, so unless the spectrum is clearly 



of the multicolour blackbody form and has no nonthermal 
tail the X-ray spectrum does not easily distinguish between 
beamed and unbeamed emission. For example, we note that 
the Galactic microquasars GRS 1915-1-105 and GRO J1655- 
40 have significantly harder emission than would be ex- 
pected from a simple multicolour blackbody (Makishima et 
al. 2000). We also point out that blazars, which are believed 
to have emission beamed towards us, tend to have relatively 
fiat uFi, spectra, varying typically by only a factor ~ 10 
from 10^" - 10^® Hz (up to an X-ray to radio ratio of 10^ for 
some Einstein Slew Survey blazars; see Fossati et al. 1998). 
In contrast, the X-ray to radio i^F,^ ratio for Galactic micro- 
quasars, which are not beamed towards us, is much greater; 
10*^ for GRS 1915-1-105 (Ogley et al. 2000 for radio; Rao et 
al. 2000 for X-ray) and 10^ for GRO J1655-40 (Hannikainen 
et al. 2000 for radio; Zhang et al. 1997 for X-ray). These 
latter sources are comparable to the brightest X-ray source 
in M82, which has an X-ray to radio ratio of at least 10^ 
(see Matsumoto et al. 2001 and Kronberg et al. 2000 for 
radio; Kaaret et al. 2001 for X-ray). Although this does not 
prove that the intermediate mass black hole candidates are 
unbeamed, it does argue in that direction. 

Additional observations would be very helpful in deter- 
mining whether beamed or unbeamed models are favoured. 
Looking at radio spectra of these sources, at high angular 
resolution, will allow a broader sample to be compared with 
blazar spectra. Another difference is the expected minimum 
time scale of variability. For an unbeamed black hole of, e.g., 
10'^ Mq, the minimum time scale would be the light cross- 
ing time across the diameter of the minimum stable orbit, 
or « 0.1 s. The time scale is likely to actually be a few times 
this; for example, Cyg X-1, with a black hole mass likely to 
be ~ 10 Mq , has a minimum timescale of significant variabil- 
ity of ~ 3 ms, three times the light crossing time (Revnivtsev 
et al. 2000). Above this timescale, the fractional variabil- 
ity increases significantly. For a beamed black hole of mass 
~ 10 Mq, the minimum time scale is decreased by a factor 
of ~ 10^ by relativistic effects and the smaller size of the 
hole itself. A long observation with Chandra in continuous 
clocking mode, with a time resolution of 3 ms, would help 
distinguish between beamed and unbeamed models. Such 
observations are important to determine whether the bright 
sources are the first representatives of a third class of black 
holes. 
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